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ABSTRACT 

The  analytical  modelling  of  croup  rosponso  phenomena  of  intervertebral  discs 
subjected  to  a constant  axial  compressive  load  is  attempted  by  using  Kolvin-solid 
modols.  A mathematical  analysis  scheme  is  proposed  for  unique  model  identification 
wherein  exact  parameter  solutions  are  developed  for  the  one-Kolvin-unit  model,  the 
threo-paramotcr-solid  model,  and  the  two-Kelvin-unit  model.  In  addition,  a method 
is  presented  by  which  the  associated  Young's  moduli  and  viscosity  coefficients  for 
an  identified  model  are  obtainable.  Hoot  importantly,  unique  parameter  values  are 
obtained  for  the  three-parameter- so lid  by  utilizing  exact  model  parameter  solutions 
on  experimental  strain,  £ (t),  data.  This  particular  model  is  observod  to  yield 
theoretical  strain,  £ (t)rn^ , values  that  aro  within  an  average  error  of  3*46$  of 
tho  experimentally  moasured  values,  6 (t)  , for  different  intervertebral  discs. 

Further,  mechanical  properties  of  the  intervertebral  discs  are  obtained  by  using 
the  values  of  the  three-parameter-solid  model  parameters  to  calculate  the  asso- 
ciated Young's  moduli  and  viscosity  coefficient.  The  appropriate  applications, 
data  limitations,  and  possible  generalizations  of  this  exact  analysis  scheme  are 


fully  discussed,  along  with  suggestions  for  future  investigatory  efforts. 


i 


V& 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH  (AF3C) 
notice  of  transmittal  to  ddc 

This  technical  report  hau  been  roviov/od  and  Is 
approved  For  public  r-l;rce  i.:d-12  (7b). 

Distribution  is  unlimited. 

A.  D.  BLOSE 

Technical  Inf  oration  Officer 


\ 


INTRODUCTION 


An  understanding  of  the  mechanical  properties  and  behavior  of  the  spinal  col- 
umn is  of  general  interest  to  researchers  in  many  areas,  ranging  from  problems  as- 
sociated with  the  man-machino  environments  in  which  man  is  subjected  to  external 
dynamic  forces  to  those  encounterod  in  the  selection  of  suitable  materials  for  in- 
tervertebral disc  replacement.  Of  particular  interest  to  this  investigation  is  the 
biomochanical  modelling  of  intervertebral  joint  response  to  constant  axial  compres- 
sive loads,  where  an  intervertebral  joint  is  understood  to  bo  an  intervertebral 
disc  and  adjoining  vertebrae. 

Nachemson  (I960),  Hircch  (1965)*  and  Rolander  (1966)  performed  experimental 
investigations  to  detormine  the  load-deflection  behavior  of  the  intervertebral  disc 
subjected  to  axial  loading.  More  recently,  Kazarian  (1975)  reported  creep  charac- 
teristics for  intervertobral  Joints  subjected  to  a constant  axial  stress,  and 
Kazarian  and  Kaleps  (1979)  illustrated  the  determination  of  Young's  moduli  and  a 
viscosity  coefficient  for  a three-parameter- solid  model  based  on  this  data. 

A basic  problem  associated  with  the  modelling  of  experimental  compressive  creep 
data  is  the  identification  of  the  most  physically  appropriate  and  analytically 
uniquo  model.  A lumped-parameter  model  consisting  of  a number  of  Kelvin  solids  in 
series  may  bo  sufficient  in  simulating  experimentally  observed  creep  phenomena. 

However,  the  number  of  Kelvin  solids  represented  in  such  a model  is  generally  de- 
pendent on  the  analytical  and  physical  intuition  of  the  investigator  and  his  wil- 
lingness and  capability  of  identifying  the  values  of  unknown,  independent  model 
parameters  via  optimization  techniques  using  the  available  experimental  data. 

One  objective  of  this  continuing  research  is  to  replace  investigative  intu- 
ition with  analytical  certitude  in  model  identification.  The  exact  mathematical 
solutions  to  a logical  progression  of  mochanical  modols  that  are  appropriate  for 
experimental  creep  data  analysis  are  presented  and  fully  discussed.  Due  to  the 
inhorent  mathematical  complexity  of  the  sxact  solution  approach  to  model  identi- 
fication for  a Kolvin  chain  consisting  of  more  than  three  Kelvin  solids,  approx- 
imating solutions  appear  warranted  for  the  aeneralizod  basic  problem.  Thn 
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presented  herein  for  throe  different  analytical  models  can  bo  easily  programed  on 
a digital  computer  to  test  experimentally  measurod  creep  data  for  unique  model 
identification  and  to  calculate  the  associated  Young's  moduli  and  viscosity  coef- 
ficients for  an  identified  model. 

ON  MODELLING  COMPRESSIVE  CREEP  PHENOMENA  WITH  KELVIN  SOLIDS 
Typical  experimental  compressive  creep  data  of  strain  £(t),  as  illustrated 
by  the  dotted  curve  in  Fig.  1 after  Kazarian  (1975),  exhibits  behavior  which  is 
similar  to  that  of  a Kelvin  solid,  where  a one-Kelvin-unit  solid  is  understood  to 
be  a spring  and  dashpot  in  parallel.  Indeed,  it  is  being  suggested  that  a Kelvin 
chain,  consisting  of  two  or  more  Kelvin  elements  connected  in  series,  night  bs 
sufficient  to  simulate  compressivo  creep  phenomena.  The  problem,  of  course,  is  to 
identify  the  minimum  number  of  Kelvin  units  in  the  chain  and  to  determine  the  in- 
dependent parameter  valuo3  associated  with  oach  Kelvin  elemsnt,  such  that  the  pre- 
dicted strain  values  of  the  Kelvin  chain  model  are  in  good  agreement  with  avail- 
able experimental  croep  data. 

The  general  solution  for  the  strain  behavior  of  an  N-Kelvin-unit  solid  is 
given  by  Fluggo  (1975)  a 3 

G(t)  -21  A.(l  - e_Xit),  (1) 

i*l  1 

which  would  appear  to  be  a logical  test  of  the  experimental  compressive  creep  data. 
When  a computerized  optimization  scheme  eiqploying  Eq.  1 was  applied  to  Kazarian's 
creep  data  for  the  T4  - T5  intervertobral  Joint,  one  value  of  H yielding  a "rea- 
sonably" good  fit  occured  for  N - 5«  A plot  of  the  calculated  strain  £ (t) 

cal 

value 3 ucing  a five-Kelvin-unit  model  is  illustrated  in  Fig.  1 by  a solid  line. 
Clearly,  after  the  first  few  minutes,  the  experimentally  observed  strain  response 
is  rather  well  predicted  by  the  fivo-Kelvin-unit  model.  Although  this  result  ap- 
pears encouraging,  further  investigation  revealed  that  the  values  used  for  the  ten 
parameters  ( 7-^,  A^;  i - 1,  2,  3,  4,  5)  were  not  unique.  That  is,  these  values 
could  be  altered  slightly,  or  even  drastically,  with  a resulting  "reasonably  good 
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fit"  to  the  experimental  creep  response  data.  This  alarming  observation  suggests 
that  another  model,  perhaps  simpler,  might  yield  a unique  fit  to  the  experimental 
data.  Onco  this  modelling  uncertainty  was  roalized,  no  further  effort  was  made  to 
improve  the  optimization  scheme  usod  to  identify  the  parameter  values  with  respect 
to  the  experimental  data.  Further,  it  is  suggested  that  the  use  of  sueh  a model, 
although  instructive,  is  of  questionable  value.  It  is  generally  felt  that  "gues- 
sing" at  parameter  values,  which  is  inherent  in  any  computerized  optimization 
scheme,  should  be  avoided,  if  at  all  possible. 

From  this  preliminary  modelling  attempt,  the  assumption  that  a Kelvin  chain 
might  be  sufficient  in  simulating  compressive  creep  response  of  intervertebral 
discs  appeared  reasonable.  The  immediate  task  was  to  uniquely  identify  the  appro- 
priate number  of  Kelvin  elements  in  the  Kelvin  chain  and  the  model  parameter  values 
by  a more  exact  mathematical  method.  To  this  end,  an  analytical  scheme  is  proposed 
wherein  oxact  parameter  solutions  are  derived  for  a few  illustrative  models.  To  be 
more  specific,  exact  parameter  solutions  aro  developed  for  the  ono-Kolvin-unit 
model,  the  threo-paraneter- solid  model,  and  the  two-Kelvin-unit  model,  along  with  a 
discussion  of  their  appropriate  applications  and  limitations  when  used  to  analyze 
experimental  strain  data. 

EXACT  PARAMETER  SOLUTIONS  TO  THE  ONE-KELVIN-UNIT  MODEL 
The  one-Kelvin-unit  model,  schematically  represented  in  Fig.  2,  is  a linear 
model,  since  the  strain  behavior  is  directly  proportional  to  the  applied  stress. 
Taking  E^  and  }{ ^ as  representing  Young's  modulus  and  the  coefficient  of  viscosity, 
respectively,  the  total  stress  applied,  0~ , can  be  expressed  by 

<r  = Exe  + n,£  * q0e  - qx€  , (2) 

where  the  dot  denotes  time-differentiation  and  the  second  equality  is  just  the 
standard  form  for  the  total  stress.  and  are  related  to  the  unit  elasticity. 


, and  damping  coefficient,  c^,  by  the  equations 


it 

y? 

v=> 

> 

(3) 

and 

Yly  S Cj  i/A, 

(4) 

*ticro  X is  the  initial  length  and  A io  the  offectivo  crosB-sectional  area  of  the 


specimen.  For  a constant  otreos 


<r  = <rQ,  t >o,  (5) 

tho  oolution  for  tho  strain,  £ (t),  to  the  constitutive  equation,  Eq.  2,  under  the 
condition  of  creep  may  bo  expressed  by 

6(t)  - (£J(t),  (6) 

where  J(t)  is  tho  creep  compliance.  According  to  Fliigge  (1975),  the  creep  coc^li- 
ance  for  this  modol  is  given  by 

j(t)  s(V<V,)(l  - e"^),  (7) 

whGre  the  parameter  is 

= qo/qi*  (*> 

Combining  Eqs.  6 and  7,  tho  solution  for  the  strain  behavior  of  the  Kelvin  solid 

can  be  represented  in  the  general  fora 

G(t)  - Al(1  - e"  (9) 

where  ^ and  A^  are  independent  parameters  requiring  determination  for  the  com- 
plete specification  of  G (t).  At  this  point,  using  Eqs.  3,  6,  7,  8,  and  9,  it  is 
easily  3hown  that 

H = VAi  do) 

and  = <rj%  i*i  • (U) 

Clearly,  Eqs.  10  and  11  can  be  utilizod  in  determing  Young's  modulus,  E^,  and  the 
viscosity  coefficient,  7^,  once  unique  values  for  the  parameters  ^ and  A^  are 
obtained. 

Assuming  this  simple  model  can  be  used  to  analyze  observed  comp re ssive  creep 
phenomena,  then  the  experimental  data  for  G{ t)  and  t may  be  expressed  in  the  form 


(12.1) 


c V1  “ ^i), 

6i*ls  VI  " '"^i*!).  (12.2) 

Even  though  these  two  equations  involve  only  the  two  unknown  parameters  and 
Ap  their  exact  solution  is  not  trivial,  unless  we  use  an  "interpolation  trick 


-ij  *■  '•  • , 
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That  is,  interpolation  of  the  physical  data  can  always  be  accomplished  such  that 


titl  ~ ti~  ti+2  “ *1*1  3 * * ' - tivn  ” tiy(n-l)  5 ^ ( 
where  y is  a "small  tine  increment".  Using  tho  requirement  of  Eq.  13  with  Eqs. 
12.1  and  12.2  and  generalizing  yields 


e.  c . *■  ^1^1/ ■>  - "i*  . — 

^1*2  “ “ A1°  1 Hi  - e 1 )«  1 » 


(14.1) 

(14.2) 


where  the  left-hand- side  of  each  equation  represents  the  difference  botween  two  con- 
secutive experimental  strain  values.  The  exact  solutions  for  the  model  parameters 
and  A^  are  inaediatoly  obtained  in  terms  of  the  experimental  strain  data  by  si- 
multaneously solving  Eqs.  14.1  and  14.2: 

= d/r)in((  - 6 l)A  eU2  - eitl»,  (15) 

A1  “ ( ^i^i+2  “ €i+l^  ^i  " 2^U1  + <16) 

It  should  bo  emphasized  that  the  model  parameter  solutions  of  Eqs.  15  and  16  arc 
dependent  on  only  the  experimentally  obtained  strain,  6(t),  data  and  can  be  easily 
esployed  by  a simplo  computer  program  to  analyze  such  data.  If  indeed  this  modal  is 
appropriate  for  the  analysis  of  any  type  of  experimental  strain  data,  then  Eqs.  10 
and  11  can  be  utilized  to  predict  the  mechanical  properties.  Unfortunately,  when 
this  simplo  modol  was  used  with  the  experimental  strain  data  of  Kazarian  (1975),  no 
unique  values  for  and  A^  were  realized.  A minimization  scheme  using  the  exact 
parameter  solutions  of  Eqs.  15  and  16  could  be  developed  for  a more  complete  analysis 
of  the  experimental  data;  however,  the  characteristics  of  the  Kazarian  (1975)  creep 
reeponco  data  suggests  that  a more  complex  model  is  required.  Since  tho  experimental 
data  graphed  in  Pig.  1 reveals  an  elastic  response  initially,  followed  by  a visco- 
elastic response,  then  perhaps  a spring  placed  in  serlos  with  a Kelvin  unit  would  be 
appropriate  as  the  next  model  to  investigate. 

EXACT  SOLUTIONS  TO  THE  THREE-PARAMETER-SQLID  MODEL 
A three-parameter-oolid  model  is  illustrated  in  Pig.  3 as  a Kelvin  unit  and 
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spring  connected  in  eorics.  Allowing  that  E^,  and  E 2 are  related  to  c^,  k^, 

and  k2,  respectively,  by  equations  similar  to  Eqo.  3 and  4,  then  the  constitutive 
equation  for  the  total  stress  applied  con  bo  derived  in  the  fora 

f + e2))^  =(e1e2/(e1  * b2))€  ♦ ( ♦ e2))6,  (17) 

where  the  dot  denotes  time-differentiation.  Eq.  17  can  be  written  in  the  normalized 
fora  as  given  by  Fluggo  (1975)* 

<r  + =%£•*■  » (18) 

by  simply  defining 

PX  5 ♦ E2),  (19) 

qoSE1E2/(E1lE2),  (20) 

qis>1iV(Ei  *V-  (21) 

The  general  solution  for  the  strain  of  Eq.  18  is  again  represented  by  Eq.  6,  only 
now  the  creep  compliance,  as  given  by  Fliigge  (1975),  is 

J(t)  = (l/qQ)(l  - e”**)  t (p1/q1)o_;tt,  (22) 

with  X-%/0^  (23) 

From  Eqs.  6 and  22,  a general  form  for  the  strain  solution  can  be  expressed  by 

£ (t)  = A^(l  - e 1 ) + k2o  ^1  , (24) 

where  the  model  paramotore  A^,  and  A2  are  given  by 

*l'qo/ql»  (25) 

*1  * ro/qG,  (26) 

(27) 

Using  these  throe  equations  with  Eqs.  19,  20,  and  21  allows  for  the  determination 
of  Young's  moduli,  E^  and  E2»  and  the  viscosity  coefficient,  Ap  in  terms  of  the 
unknown  model  parameters: 

El  = ^/(Ax  - A2),  (28) 

E2=VA  2»  (29) 


ni  * ^/\(A1  “ A2>  = V^l* 


As  before,  valueB  for  the  model  parameters  must  bo  uniquely  determined  by  analyzing 
the  experimentally  available  strain  data. 

Uoing  Eq.  24  as  the  model  "data  equation",  then  any  threo  consecutive  strain 
data  values  obeying  the  interpolation  requirement  of  Eq.  13  can  be  represented  by 


6i  - Ax(l  - * A2e";ilti, 

€i+l  * Ax(l  - o 1 ie  1 ) ♦ A2e  1 *-0  1 , 

C _ . /1  - X^t^  -2  . - ^l^i  *2 

ti#.2  SV1  - e 0 A ) ■*•  A2e  A xe  1 . 


(31.1) 


(31.2) 


(31.3) 


Exact  solutions  for  the  three  model  parameters  are  now  obtainable  in  terms  of  the  ex- 
perimental strain  data  by  solving  the  above  throe-equations  simultaneously.  The 
results  for  2-^  and  A^  are  identical  to  those  given  by  Eqs.  15  and  16,  respectively, 
and  A2  can  be  expressed  by 

A2  = - AjU  - o'  Vi)).  *1*1.  (32) 

Kazarian  and  Kaleps  (1979)  recently  analyzed  compressive  creep  phenomena  using 
a three-paramo tor- solid  modol  described  by  an  equation  like  Eq.  24.  Their  results 
for  the  Young' b moduli  and  the  coefficient  of  viscosity  for  a few  intervertebral 
Joints  aro  given  in  Table  1 as  E*,  E*,  and  By  their  analysis  scheme,  which  em- 

ployed various  approximations  and  optimization  schemes,  the  predicted  strain  values, 

£ (t  ) , compared  favorably  with  the  experimental  strain  data,  £ (t.)  , ranging 

i cal  l exp 

from  about  1%  to  15$  with  about  a 5$  average. 

When  the  exact  parameter  solutions  developed  herein  were  employed  to  analyze 
the  Kazarian  (1975)  data,  unique  parameter  valuos  were  obtained  for  the  interverte- 
bral Joints  considered  in  Table  1.  The  results  for  the  Young's  moduli  and  the  vis- 
cosity coefficient  are  indicated  in  Table  1 as  E^,  E^,  and  and  agree  favorably 

with  the  results  predicted  by  Kazarian  and  Kaleps  (1979).  When  the  first  three  ex- 
perimental strain  values  (elastic  reeponse)  are  eliminated  from  the  data  points  as- 
sociated with  each  spinal  cegnent,  the  exact  analysis  scheme  prodicts  strain  values. 


6 (t, ) that  arc  within  an  average  of  3.48  % of  the  experimental  data,  £ (t. ) , 

i.  c&j.  i exp 

for  the  twenty  epinal  segments  considered.  For  each  spinal  segment  considered,  the 
"average  of  the  absoluto  % error,  £,"  is  defined  by 

1 1 u/»> Jb  ABsfu  - €(\)ci (33> 

and  is  observed  to  bo  rather  small  with  the  exception  of  test  I.  D.  No.  1 and  16. 
With  these  two  Epinal  segments  eliminated,  the  predicted  strain  values,  £ ^i^cal* 
are  within  an  average  error  of  2.21  % of  the  experimental  data,  £ (t. ) , for  the 

X €Xp 

remaining  eighteen  spinal  segments.  For  these  intervertebral  Joints,  a graph  of 
strain  versus  tina  for  £ (t^cai  will  b®,  essentially,  superimposed  on  the  corres- 
ponding graph  of  £ (tj^exp*  excePt  for,  approximately,  the  first  threo  strain  values. 

Inclusion  of  tho  first  threo  predicted  strain  values  for  each  spinal  segment 
results  in  an  average  discrepancy  of  8%  between  £ (t. ) 1 and  £ (t. ) for  the 

twenty  intervertebral  Joints  of  Tablo  1.  It  should  be  noted  that  tho  exact  param- 
eter solutions  are  very  sensitive  to  the  experimental  strain  values  and,  thus, 
critically  dependent  on  the  smoothing  and  interpolation  programs  employod  to  reduce 
the  original  data.  Certainly,  an  optimization  scheme  utilizing  calculated  param- 
eter values  could  bo  developed  for  a more  complete  analysis  of  the  experimental 
data,  in  an  attempt  to  obtain  closer  agreement  between  £ (t. ) . and  £ (t  ) for 

X Cal  X 6Xp 

the  first  three  data  points  of  each  intervertebral  joint.  However,  the  compressive 
creep  datum  analyzed  tends  to  suggest  that  a slightly  more  complex  model  may  be 
required  for  improved  agreement  between  theoretical  and  experimental  predictions. 


EXACT  PARAMETER  SOLUTIONS  TO  THE  TVO-KELVIN-UNIT  MODEL 
The  constitutive  equation  for  tho  two-Kolvin-unit  solid,  as  illustr&tod  in 


Pig.  U,  can  be  expressed  in  the  normalized  form  ae 

<r+  p =q0£  f qx£  + q2d  » (34 

where  the  dot  denotes  time-differentiation.  The  stress  and  strain  coefficients  of 
Eq.  34  are  related  to  Young's  moduli,  E^  and  E^ , and  the  viscoeity  coefficients, 
and  ??2»  b7  tho  following  equations; 


P ! r ( +■  V2)AE1  + E2)»  (35) 

% = ElV(El  + V'  (36) 

qx  = (^!E2  + n2\)/{£l  * E2)#  (37) 

q2  = ^l^V^l  * V*  (36) 

The  general  {solution  to  Eq.  34  for  the  strain,  as  given  by  Fliiggct  (1975),  is 

6 (t)  = ^J(t)  = ro(l  + Pl^)(l  - •■*lt)/q2*i(*2  - \) 

* <roU  + P1^2)(l  - «'^)/q2VH  - ^),  (39) 

where  and  ^ aro  the  roots  of 

q2*2  “ qlX  + % ' 0>  (4°) 

Taking  as  the  negative  root  and  /t-,  as  the  positive  root  of  Eq.  40,  and 

4-^  5 (^(^  h Pl^l^q2^l^2  ""  ^2.) * (41) 

A2  c ro(1  * P1^2^q2^2(:Ll  “ (42) 

then  the  general  solution  for  the  strain  behavior  can  bo  expressed  by 

6 (t)  = A^(l  - e 1 ) + A2(l  - e ^2  ),  (43) 

Young's  moduli  and  the  viscosity  coefficients  can  be  related  to  the  four  parameters 
ob  Eq.  43  by  using  Eqs.  35,  3b,  37,  33,  40,  41,  and  42: 

E1  = roa2  “ *1)/(A1(*1  ♦ V + 2^^)*  (44) 

E2  = £<*1  - + h)  * ^lV'  (45) 

^1  5 V*l*  (46) 

9*2  (47) 


Again,  value 3 for  the  four  model  parameters  must  be  uniquely  identified  by  analyzing 
the  experimentally  obtained  comprsssivo  croep  data. 

Taking  Eq.  43  as  representative  of  the  "experimental  data",  then  the  four  data 
equations  necessary  for  the  exact  solution  of  the  model  parameters  are  of  the  form 
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where  the  "tine  interval"  requirement  of  Kq.  13  has  been  incorporated.  Solving  Eqs. 
48.1  to  48.4  simultaneously  and  defining 

Ei+j  £ £ ifj  " £ i+( j-1)’  J s 1»  2»  3>  4>  (A9) 

gives  the  paramoter  solutions  in  the  form 
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0 5 (EU2  " EUlEi+3)e  + (EivA.4  - Ei*2Ei+3)o  + (Ei+3  - Ei+2Eit4)^(52) 
where  is  the  ncgativo  root  and  ^ I3  the  positive  root  of  Eq.  52. 

When  the  compressive  croep  data  of  Kazarian  (1975)  was  analyzed  by  the  above 
parameter  equations,  no  unique  set  of  values  for  Ap  Ap  and  A2  wao  obtained. 
Some  calculated  values  appeared  repetitively  as  the  experimental  data  v/as  tested, 
but  it  wa3  not  possible  to  identify  a "unique  sot"  of  parameter  values  for  any 
intorvertebral  Joint.  As  such,  use  of  the  two-Kelvin-unit  solid  as  an  analytical 
model  was  seriously  questioned. 


CONCLUSION 

The  analytical  modelling  of  comprossivo  creep  data  by  Kelvin  solids  has  been 
considered.  Unique  model  parameter  solutions  are  accomplished  by  a mathematical 
scheme,  which  has  boon  fully  illustrated  for  the  two,  throe,  and  four-parameter 
solids.  The  first  and  last  models  discussed  are  not  capable,  by  themselves,  of 
predicting  observod  compressive  creep  phenomena;  however,  successful  modelling  of 
the  experimental  data  appears  possible  by  using  the  three-parameter- solid  model  or 


moro  complex  Kelvin  chains.  Koro  precisely,  the  five-parameter- solid  model  is  sug- 
gested as  tho  next  logical  model  to  develop.  It  is  expected  to  yield  a better  fit 
to  tho  experimental  compressive  creep  data  for  the  entire  time-domain  than  the  three 
parameter- solid  model.  Exact  parameter  solutions,  of  the  type  discussed  herein, 
have  not  been  obtained  for  a Kelvin  chain  consisting  of  three  or  moro  Kelvin  element 
It  is  suggested  that  any  such  derivational  attempts  will  be  intimidated  by  the  math- 
ematical complexity  and  uncertainty  arising  from  the  algebraic  solutions.  It  should 
however,  be  possible  to  develop  a generalised  approximating  scheme,  utilizing  the 
equations  presented  herein,  which  is  capable  of  identifying  a unique  model  for  tho 
prediction  of  compressive  creep  response.  Further,  it  should  be  realized  that  ex- 
perimental data  analysed  by  this  type  of  method  should  be  carefully  smoothed  and  in- 
terpolated, because  of  the  sensitivity  of  the  equations  representing  tho  exact  so- 
lutions for  the  model  parameters. 

Tho  implications  of  this  report  are  that  with  improved  computer  programs  per- 
taining to  interpolation,  xaoothing,  and  analycic  modols  (development  of  the  five- 
par&meter-solid  modol),  all  of  Dr.  Leon  E.  Kazarian's  experimental  data  on  the  load- 
deflection  behavior  of  intervertebral  joints  can  be  quickly,  efficiently,  and  ac- 
curately reduced  via  exact  parametric  solutions,  with  associated  Young's  moduli  and 
viscosity  coefficients  obtained.  Such  an  accomplishment  would  represent  the  mo st 
authorative  ro search  results  cn  intervertebral  joints  to  date  and  should  become  the 
primary  reference  for  investigators  associated  with  problems  ranging  from  the  man- 
machino  environments  to  those  encountered  in  tho  selection  of  suitable  materials 
for  disc  replacement. 
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Figure  1 


A comparison  of  the  experimental  compressive  creep  response 
for  the  human  T4  - T5  intervertebral  Joint,  by  Kazarian,  with 
the  predictions  of  a 5-Kelvin-unit  model,  by  Burns  and  Kaleps. 


Table  1 


A comparison  of  the  Kazarian  and  Kaleps  (1979)  Young's  moduli 
(E£,  E*)  and  viscosity  coefficient  (?£*)  for  intervertebral 
joints  is  made  with  those  predicted  by  Bums  and  Kaleps  and 
represented  by  E,,  E2,  respectively.  A comparison  of  the 

predicted  strain  values,  £ (t^)  of  the  611X03  And  Kalep3 
exact  analysis  scheme  with  the  strain  values  obtained  experi- 
mentally, £ (tj.)  , by  Kazarian  (1975)  is  represented  as  an 

"average  of  tho  absolute  % error",  where  the  first  three  strain 
data  values  for  oach  intervertebral  joint  has  been  eliminated. 
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